Itinerant Ferromagnetism in an Ultracold Atom Fermi Gas 
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We address the possible occurrence of ultracold atom ferromagnetism by evaluating the free energy 
of a spin polarized Fermi gas to second order in its interaction parameter. We find that Hartree- 
Fock theory underestimates the tendency toward ferromagnetism, predict that the ferromagnetic 
transition is first order at low temperatures, and point out that the spin coherence time of gases 
prepared in a ferromagnetic state is strongly enhanced as the transition is approached. We relate 
our results to recent experiments. 
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Introduction — Itinerant ferromagnetism is common in 
metals. Nevertheless, because it flows from a strong- 
coupling Fermi liquid instability, the microscopic physics 
that controls its occurrence is less well understood than 
the physics that controls superconductivity |]|. In the 
electron gas case, for example, accurate quantum Monte 
Carlo calculations suggest |2( that the transition to the 
ferromagnetic state occurs at a critical density nearly 
three orders of magnitude smaller than that predicted by 
mean-field (Hartree-Fock) theory. Even in the simplest 
model of interacting electrons, the single-band Hubbard 
model, solid predictions on the occurrence of ferromag- 
netism are rare and often restricted to particular band 
fillings 0,0. Understanding the nature of the paramag- 
netic to ferromagnetic phase transition, when it occurs, 
has also been challenging. Experimental progress has 
recently been achieved by applying hydrostatic pressure 
to itinerant ferromagnets with a low Curie temperature, 
making it possible to study the transition in the zero 
temperature limit and test for theoretically predicted 
quantum critical 0| behavior. In these experiments, the 
line of continuous transitions in the temperature-pressure 
phase diagram appears [E H> El to terminate at a tri- 
critical point with decreasing temperature, connecting 
with a low-temperature line of first order transitions. In 
mean-field theory, first order magnetic transitions can fol- 
low from a non-monotonic quasiparticle de nsity -of-states 
near the paramagnetic state's Fermi energy [lfj. Belitz et 
al. 01 have argued, however, that coupling of the order 
parameter to gapless modes leads to nonanalytic terms in 
the free energy, and generically drives the transition first 
order. These nonanalytic terms were first pre dicted by 
Misawa on the basis of Fermi-liquid theory |l2( , and are 
a consequence of gapless particle-hole excitations. Theo- 
ries of the phase transition are still qualitative however, 
and detailed experimental corroboration of this picture 
is still lacking. 

In this Letter we address the possible complemen- 
tary realization of ferromagnetism in ultracold fermionic 
atoms, which are accurately described by a short-range 
interaction model [Til fl4j . In Hartree-Fock theory 01 



the zero-temperature ferromagnetic transition of this 
model is continuous and the ground state is ferromag- 
netic when the gas parameter, i.e., the product of the 
Fermi wave number kp of the unpolarized system, and 
the s-wave scattering length a of the short-range poten- 
tial, satisfies k-pa > w/2. The phase separation predicted 
by Houbiers et al. |l5j at the same gas parameter is one 
plausible manifestation of ferromagnetism but, as we dis- 
cuss below, not the most likely one. Trapped-atom mo- 
tivated inhomogeneous generalizations of these Hartree- 
Fock theories have recently been analyzed by Salasnich 
it et al. ^(| and Sogo and Yabu 0] . 



The issue of ferromagnetism in a two-component 
atomic Fermi gas is of particular interest because of the 
ongoing experimental study of stron gly interact ing, de- 
generate, fermionic alkali atoms [lilTlii. 120L Ell Ea E3 . 
|2 1| . The focus so far has been on observin g th e formation 
of a fermion pair condensate [2^, |2(| W\ Eil Efl l30t EH 
E3,EEE3| in the BCS-BEC crossover 0013 regime 
close to a Feshbach resonance jlM Ell ■ Our interest is in 
the repulsive interaction side of the resonance, where we 
believe it will be possible to achieve unprecedented ex- 
perimental control over ferromagnetism. In making this 
assertion we are assuming that the formation time of the 
molecular BEC state (which occurs under the same con- 
ditions when the state is prepared by crossing from the 
attractive interaction side of the resonance) can exceed 
experimental time scales when the state is prepared by 
approaching the resonance from the repulsive interaction 
side. 



The character of the ferromagnetic state that can be 
realized experimentally in these systems requires some 
comment |4(j. Since s-wave scattering does not occur 
between identical fermions, interaction effects require the 
presence of two hyperfine (pseudospin) species. Using 
standard techniques, the atomic system can be prepared 
in a pseudospin coherent (ferromagnetic) state, in which 



2 



all atoms share the same spinor: 



vac) 



k|<23fc F 



( c k a crea t es an atom with momentum k and hyperfine 
spin a.) In Eq. ip specihes the orientation of the mag- 
netic order parameter in the x — y plane and AE is the 
Zeeman energy difference between the hyperfine states. 
(Since the number of atoms in each is conserved, we can 
transform to a rotating wave picture and let AE — > 0.) 
Overall spin polarizations in the z direction are not ac- 
cessible. This fully spin coherent state always has a 
lower energy than the phase-separated state discussed 
m Refs. 

since, in the magnetic language, the 
latter has a domain wall which costs finite energy. Fer- 
romagnetism in these systems will be manifested by per- 
sistent coherence between hyperfine states. 

In this Letter we argue that ferromagnetism occurs on 
the repulsive interaction side of a Feshbach resonance. 
Our principle results are summarized in Fig.^and Fig.|3 
We find that i) Hartree-Fock theory underestimates the 
tendency towards ferromagnetism 0, u) that the tran- 
sition between ferromagnetic and paramagnetic states is 
first order at low temperatures and, in) that the coher- 
ence decay rate decreases rapidly as the thermodynamic 
stability region of the ferromagnetic state is approached 
from the repulsive side of the resonance. 
Second order perturbation theory — It is convenient to 
view the gas as a mixture of two independent noninter- 
acting gases of spinless fermions. The grand-canonical 
Hamiltonian of the system is then 

5- - M« ) ^a(x) 

={+,-} 



H= d 



g / dx^t(x)Vl(x)V-(x)V+(x) , (2) 



with g — Anati 2 /m. The chemical potentials are deter- 
mined by n a = dpoa/d/j, a , where n a is the density of 
atoms in hyperfine state \a), and the pressure of the non- 
interacting gas is given by 



POc 



V 



k 



In 1 



-P(e k -Hc.) 



(3) 



with k^T the thermal energy, V the volume, and e k = 
fi 2 k 2 /2m the single-particle dispersion. The entropy den- 
sity is determined by s = d(po + +po-)/dT, and the total 
free energy density is given by /(n+, n_) = e — Ts, with 
the total energy density expressed as the sum of three 
contributions, e = e^ + + e^ 2 \ The first two contri- 
butions correspond to Hartree-Fock theory, and are given 
by 



,(o) 
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N k „e k 
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(4) 
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FIG. 1: Magnetization £ as a function of kpa, for various tem- 
peratures. From left to right T/T F = 0,0.1,0.15,0.2,0.25. 
The dashed lines indicate magnetization jumps. The inset 
shows the critical temperature as a function of the gas pa- 
rameter. The solid line indicates first-order transitions, and 
the dotted line second-order transitions. The dashed line is 
the Hartree-Fock theory result. 



where -/V kQ , is a Fermi occupation factor. The contribu- 
tion to the energy density that is second order in inter- 
actions is given by [l3| 



,(2) 



V^/iV kl)+ iV k2 ,_(iV k3 , 

1/3 2^ 



£ki + Ck 2 



£k 3 



£ k 4 



(5) 



where the prime indicates that the sum is over wave vec- 
tors such that ki +k 2 = k 3 +k4. The above second order 
correction takes into account the so-called unitarity limit, 
i.e., the energy dependence of the vacuum scattering am- 
plitude to all orders in ka, to second order. 
Results — The magnetization results, summarized in 
Fig.Q were obtained by numerically minimizing the total 
free energy /(n_,n+) vs. £ = (n+ — n_)/(n+ +n_), for 
a series of temperatures and total densities n + + n_ = 
k F /3ir 2 . At zero temperature, we find that the system be- 
comes partially polarized if kpa > 1.054, and reaches the 
fully-polarized state at kpa — 1.112. For higher tempera- 
tures interactions have to be stronger to polarize the sys- 
tem. For temperatures T < T tc , where T tc ~ 0.2Tp with 
Tp the Fermi temperature, the transition is discontinu- 
ous, and the magnetization exhibits a jump. The jump 
becomes smaller with increasing temperature, vanishing 
at Tt c . The inset shows the transition temperature as a 
function of kpa. A line of first-order transitions, denoted 
by the solid line, joins a line of continuous transitions, 
denoted by the dotted line at T = T tc and kpa — 1.119. 

The first order behavior at low temperatures is ex- 
pected on the basis of the arguments of Belitz et al. • 
In our case the gapless modes that drive the transition 
first order are particle-hole excitations. The coupling 
of these excitations to the magnetization is neglected in 
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Hartree-Fock theory, which therefore always predicts a 
continuous transition. Eq. JSJ takes the coupling between 
the magnetization and the particle-hole excitations into 
account to lowest order. 

Experimental implications — The ferromagnetic state 
can be identified by measuring the interaction energy, ei- 
ther by studying the expansion properties of the gas [41| , 
or by using RF spectroscopy |Hi S E3 ■ The full y 
polarized state is distinguished by the absence of any in- 
teraction energy. In the experiments by Bourdel et al. 

on 6 Li gases, the interaction energy appears to van- 
ish when the regime of strong repulsive interactions is ap- 
proached. In Fig.|21we plot the interaction energy divided 
by the kinetic energy for their experimental parameters, 
as a function of the magnetic field. Given the fact that 
we have not taken into account the inhomogeneity of the 
system, the agreement is remarkable, strongly suggesting 
that a ferromagnetic transition occurs in this system. If 
we interpret the experimental data accordingly, the tran- 
sition is found to occur at fcpa ~ 1 at T = 0.6Tp, which is 
slightly smaller than our calculated value (kpa = 1.56) at 
this temperature. In the experiments of Gupta et al. |33| 
the value of fcpa at which the mean-field shift appears to 
vanish is even smaller compared to the value we predict 
for the onset of ferromagnetism. Since the atom system 
in these experiments is prepared in a ferromagnetic state, 
these discrepancies between theory and experiment could 
be due to the rapid increase in spin coherence time which 
is expected as stable ferromagnetism is approached, as we 
now explain. 

Pseudospin decoherence in these systems is due to spa- 
tial inhomogeneities in the Zeeman energy AE. Suppose 
the potential the | f) atoms feel is E^(x), and the poten- 
tial the | |) atoms feel is E^(x). The decay rate of the 
fully coherent state is suppressed because the quasiparti- 
cle energies of the unoccupied pseudospins are shifted by 
the interactions. Fermi's golden rule implies a coherence 
decay rate 

r = y^ I A£ k ,, k | 2 <S(e k - e k , - gn) (JV k ,+ - N Vi _) , 

k'.k 

(6) 

where 

A^ k = l|^( Et ""; E 'W )e^-»-. (7) 

The consequences of interactions can be illustrated by 
taking |Al? k / ik | 2 = 5E 2 e~ A ( k ~ k ') where A is the length 
scale of magnetic-field inhomogeneities. In Fig. [3] the 
T = spin coherence time is shown for a series values 
of A. These results were obtained by taking a trapping 
frequency oj/2tt = 20 Hz , estimating SE as the difference 
in Zeeman splitting change between the edge and center 
of the cloud, and assuming 7 x 10 6 atoms at a density 
n~4x 10 13 cm -3 , following Gupta et al. jl^- Clearly, 
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FIG. 2: Interaction energy divided by kinetic energy as a 
function of magnetic field for the experimental parameters 
of Bourdel et al. |4l|. We take a temperature T = 3.5 /xK 
= 0.6Tf. For details on the magnetic-field dependence of the 
scattering length see, for example, Ref. [4ll |. 
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FIG. 3: Magnetic field inhomogeneity limit on the spin co- 
herence time of the fully-polarized state. 

the spin coherence time is strongly enhanced for increas- 
ing interactions. The difference in experimental results 
between Bourdel et al. ^lj an( A Gupta et al. .'i.'i might 
be related to differences in magnetic field. We note that 
the magnetic-field inhomogeneities are necessary for the 
equilibration of the hyperfine spin degrees of freedom. 
Since molecule formation cannot occur in the fully po- 
larized state, ferromagnetism competes kinetically with 
Bose-Einstein condensation of molecules 0] • 

Coherence decay and atomic ferromagnetism can also 
be studied by measuring the size of the cloud. In a local 
density approximation, valid since the oscillator length 
exceeds the Fermi wavelength, the size of the cloud is 
proportional to the square root of the Fermi energy. It 
follows that the radius of the fully-polarized state is a 
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factor 2 1 / 3 larger than that of the unpolarized state. 
The first order character of the phase transition could be 
detected by performing experiments with a mixture of 
fcrmions and bosons. (The interactions between bosons 
and the fermions should be weak enough to make boson 
mediated attractive interactions between fcrmions negli- 
gible.) Suppose for example that the mixed system is in 
equilibrium, and that the energy and number of atoms 
are conserved as the bias field is varied. Adiabatically in- 
creasing k-pa from the paramagnetic to the ferromagnetic 
regime will lead to a temperature increase that is tied to 
the entropy reduction in the ordered state. For T < T tc 
the temperature variation should be hysteretic. These 
temperature changes, although typically relatively small 
(~ 10~ 3 Tp), are larger for a smaller boson to fermion 
mass ratio and boson concentration and might be ob- 
servable. 

Discussion and conclusions. - Although ferromag- 
netism is a strong coupling instability and our theory 
is perturbative, we nevertheless believe that the phase- 
diagram in the inset of Fig.^is reliable. The interaction 
energy of the fully-polarized state, which is an eigenstate 
of the full hamiltonian, is exactly zero. Moreover, a cal- 
culation to third order in the gas parameter shows that 
the energy of the paramagnetic state energy is increased 
in comparison to the second-order result |43| . Hence we 
expect that the second-order perturbation theory under- 
estimates the transition gas parameter. However, since 
consistency requires that the critical point lies in the 
strong-coupling regime where k-pa ~ 1, there appears 
to be little room for movement. An experimental de- 
termination of the phase diagram appears to be within 
reach and would be interesting. The magnetic properties 
of ultracold Fermi gases could provide a very interest- 
ing system to explore fundamental aspects of ferromag- 
netism, including the dynamics of domains walls which 
could be directly manipulated by a one-way barriers , 
the noncquilibrium formation of the ferromagnetic state, 
spin waves, and spin transfer effects. 
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